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Summary 


The North American drought in the summer of 1988 was one of the worst 
droughts of the last 50 years. The rainfall was 50 to 85 percent below normal in 
areas of the Midwest, northern Plains and Rockies. This drought caused an 
estimated $40 billion in damages in the United States, as a result of damaged 
crops, reduced river flow, widespread forest fires and record high surface 
temperatures. The primary cause of this drought has previously been identified as 
abnormally cool ocean temperatures over the eastern Pacific Ocean, called a La 
Nina” (the opposite counterpart to the well-known El Nino). However, despite 
the use of observed ocean temperatures from 1988, this drought has not been well 
reproduced by any global-scale climate model. The GEOS model at GSFC did 
simulate some of the features of a drought over North America, but the 
accompanying circulation was not so well simulated, leading one to infer that the 
model had simulated the drought for the wrong reason. To isolate the causes for 
this remarkable circulation, wind, soil moisture, and vegetation data, taken from a 
blend of satellite observations and data analysis, were used to steer the model to 
yield a better prediction. This study found that many of the models deficiencies 
result from wind biases that are carried into North America from surroundings 
regions. This is not an unexpected result because weather and climate have global 
connectivity, and model errors tend to propagate from one region to another. 
Through such errors in the wind, the transports of heat and moisture into the 
North American region were changed. This study also found that the model is 
improved more through the use of soil moisture data than by vegetation data. In 
this way, we have understood what aspects of the simulation lead to what type of 
prediction errors. 
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On Simulating the mid-western-US Drought of 1988 with GEOS-2 GCM 


ABSTRACT 

The primary cause of the midwestem North American drought in the summer of 1988 has been 
identified to be the La Nina SST anomalies. Yet with the SST anomalies prescribed, this drought 
has not been simulated satisfactorily by any general circulation model. Seven simulation- 
experiments, each containing an ensemble of 4-sets of simulations, were conducted with the 
GEOS GCM for both 1987 and 1988. All simulations started from January 1 and continued 
through the end of August. In the first baseline case, Case 1, only the SST anomalies anc some 
vegetation parameters were prescribed, while everything else (such as soil moisture, snow-cover, 
and clouds) was interactive. The GCM did produce some of the circulation features of a drought 
over North America, but they could only be identified on the planetary scales. The 1988 minus 
1987 precipitation differences show that the GCM was successful in simulating reduced 
precipitation in the mid-west, but the accompanying circulation anomalies were not well 
simulated, leading one to infer that the GCM has simulated the drought for the wrong reason. To 
isolate the causes for this unremarkable circulation, analyzed winds and soil moisture were 
prescribed in Case 2 and Case 3 as continuous updates by direct replacement of the GCM- 
predicted fields. These cases show that a large number of simulation biases emanate from wind 
biases that are carried into the North American region from surroundings regions. Inclusion of 
soil moisture also helps to ameliorate the strong feedback, perhaps even stronger than that of the 
real atmosphere, between soil moisture and precipitation. Case 2 simulated one type of surface 
temperature anomaly pattern, whereas Case 3 with the prescribed soil moisture produced another. 
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In the combined case. Case 4, the two patterns merged and produced somewhat more realistic 
evapotranspiration and precipitation patterns. 

Case 5 had no soil moisture and winds updates inside a North American region, but with these 
updates outside this region. The simulation in Case 5 is very similar to Case C4, and shows that 
many of the local simulation errors originate outside the North American region. This is not an 
unexpected result because if weather and climate have global connectivity, errors will propagate; 
however, such a large connectivity at a seasonal scale is a new result. One expects it would not 
be too different for other state-of-the-art climate models. 

Cases 6 and 7 were built on Case 5, with the addition of correct and opposite-year vegetation data 
from ISLSCP, respectively. The results show that the model is unable to reveal discernible 
improvement in the correct vegetation simulations for the drought year. One can conclude that 
the soil moisture governs the outcome of the land-atmosphere feedback interaction much more 
than the vegetation parameters. Clearly, vegetation and soil moisture go together, but if the soil 
moisture effect dominates, then its biases would mask the influence of corrected vegetation. 

The authors argue that model biases strongly influence the prediction errors. Through wind- 
errors, they change the transports of heat and moisture into the North American region. When 
winds were prescribed from analysis of observations outside the chosen region, the model 
produced a much better JJA drought as compared to the control, and it remained almost as good 
as the one in which winds and soil moisture were prescribed everywhere. This shows that biases 
in circulation and advective transports propagate and strongly contribute to the simulation biases. 


1. Introduction 
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The Goddard Laboratory for Atmospheres has a well-documented state-of-the-art general 
circulation model (GCM) called the GEOS GCM (Conaty et al., 2001) that is usually run with 
coupled land and prescribed sea-surface temperatures. A climate version with only 20-sigma 
layers in the vertical and several physics upgrades is used at the Climate and Radiation Branch 
(CRB) of the Laboratory for Atmospheres. Specifically, the CRB version of the GEOS GCM 
employs a Simplified version of the Simple Biosphere Model (SSiB, Xue et al., 1991) for land 
processes that is further modified with snow physics package developed by Mocko and Sud 
(2001). The cloud physics consisting of Microphysics of Clouds with Relaxed Arakawa- 
Schubert (McRAS) due to Sud and Walker (1999a&b) handles all moist processes. Such a 
model simulates climate, climate change, and its biogeophysical consequences consistently even 
if it does not capture all of the features of a specific climatic episode under investigation. Thus, 
in interpreting a model’s forecast, its realism is often an issue whereas its internal consistency is 
not. There can be several causes for a model’s failure to accurately simulate an observed climatic 
episode. Among them are coarse resolution, simplifications in the representation of atmospheric 
physics (the primary cause of intrinsic model deficiencies particularly due to parameterizations), 
and natural variability of the actual as well as observed climate system. The natural variability of 
climate is a major source of unpredictability. Consequently, one must view observations as a 
single realization amongst a host of possible climate pathways in nature that a particular initial 
state might have produced. Regardless, as research tools, Earth-system models can help us 
understand and discern the roles of coupled land-atmosphere-ocean interactions in maintaining 
and modulating the evolving climate of the Earth, which includes major hydrological events such 
as droughts and floods. 

The drought of the summer of 1988 over the mid-western United States was a major North 
American drought. This drought persisted over the agricultural region of the Great Plains during 
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the spring and early summer, and had a devastating effect on crop-yields in the mid-west as well 
as the US economy as a whole (Trenberth and Branstator, 1992). Its catastrophic features 
included: i) 50-85% below normal precipitation in the midwestem North America, Northern 
Plains, and Rockies; ii) record high surface temperatures; widespread forest fires that burnt 
nearly 4.1 million acres of forests by mid-autumn; iii) record low Mississippi river discharge 
(40% of normal in mid-June 1988); and iv) a total estimated economic loss of roughly 40 billion 
dollars. On the basis of this data, the drought of 1988 turned out to be the worst drought in the 
last 40 years. It was accompanied by surface flux-anomalies on the order of about 60 W m' , and 
yet it is this drought that the GCMs of our times often fail to simulate. Land feedbacks such as 
surface albedo (Sud and Molod, 1988), soil moisture (Dirmeyer 1999), and vegetation variations 
(Sud et al., 1993, 1995) have been shown to produce a positive feedback on a drought circulation 
and rainfall; however, they do not help much in explaining or understanding transient droughts. 
Some successes with regional models ran with observed lateral data have been documented (e.g., 
Giorgi et al., 1996; Hong and Pan, 2000; Jenkins and Barron, 2000), but the question remains: 
how does one obtain reliable lateral forcing data to predict such a drought? 

With every major model improvement of the GEOS GCM (and its earlier versions), we 
attempted to simulate the North American summer drought of 1988, but thus far have had only 
limited success (e.g., Mocko et al., 1999). Naturally, these failures have provided a daunting 
challenge to really determine whether it is the model, or boundary forcings like soil moisture, or 
the hard-to-simulate pathway that the climate system chose. Is the soil moisture, or the model’s 
biases, or remote forcing errors, or poor representation of convection the key contributing factor? 
After nearly a 14-year time-gap, we are still expecting to find a good explanation for the colossal 
lack of predictability of this drought. It would be interesting to determine what is missing in the 
model and how various aspects of model deficiencies might interact with each other to cause the 
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model to not simulate the drought. In this study, we shall explore the key observational features 
of the drought of 1988 and the primary reasons for lack of (or limited) success in simulating it. 

The Climate Analysis Center at National Center for Environmental Prediction (NCEP, formerly, 
National Meteorological Center, NMC) produced the following schematic picture (Fig. 1) to 
hi ghligh t some key dynamic features of the drought of 1988. It shows how the drought 
perpetuated. The jet stream in the mid-west was located far north of its normal position in 
association with an anomalous ridge of high pressure in the northern plains. This led to 
northward transport of moist air masses along the West Coast of North America making the 
region rainy and damp. In addition, the high-pressure system over the Great Plains caused the 
low-level jet (normally bringing the Gulf moisture into the Great Plains region) to weaken, 
thereby shutting off the moisture supply to the region. In fact, this low-level jet is so important 
that in all the other recent droughts over the mid-west North America, e.g., summer of 2000, 
weakening of the low-level jet producing a reduction in moisture transport was a key factor in the 
production/maintenance of the droughts. The accompanying sea-level pressure picture for June- 
July- August (JJA) of 1988 minus 1987 (Fig. 2) shows high pressure over Canada was primarily 
responsible for bringing in dry air from the north. Together with a cyclone over the Northern 
Rockies and an anticyclone to its south, one notes a weakening of the low-level jet. In this way, 
the low-level jet carrying moisture laden low-level Gulf air was replaced with dry air from 
Northern Canada. This scenario would naturally lead to a drought. The drought of 1988 
persisted through the early summer and then started to fade away in the late summer when the 
low-level jet strengthened again and copious rains returned to the region. Further analysis of the 
physical processes of the drought can be found in Trenberth and Guillemot (1996). The 
normalized difference vegetation index (NDVI) analysis (Fig. 3a) identifies the drought region 
along the low-level circulation of dry airmass transports from Canada. This circulation produces 
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moisture divergence and drying in the drought region. The drought had its greatest impact in the 
northern Great Plains. It intensified over the northern Great Plains and spread across much of the 
eastern half of the United States with total precipitation for April through June of 1988 being 
lower than the Dust Bowl period. The summer of 1988 is well known for the extensive forest 
fires that burned across western North America, including the catastrophic Yellowstone fire. In 
addition to dry conditions, heat waves during the summer of 1988 broke long-standing 
temperature records in many mid- western and northeastern metropolitan areas. 

Because the mid-western drought of summer 1988 was associated with a cold SST Event over 
the Tropical Oceans, popularly called La Nina, it is natural to explore its link to the drought. A 
40 years dataset (1947-1986) from COAPS analysis (Sittel, 1994) was examined. It has 
identified Great Plains regions to be susceptible to warming/droughts in association with cold La 
Nina episodes (Figure 4). Castro et al. (2001) examined the NCEP-NCAR reanalysis (1948- 
1998) and showed that La Nina conditions tend to shorten the spring season rainfall in the Great 
Plains and lead to drier early-summer conditions. Such evaluations of climatic variables of the 
Northern Hemisphere as a whole discern the sequence of phases during the ENSO cycles. Some 
applications of this investigation could be used in preparation for ENSO-related drought and 
forest fires as well as provide differential guidance for GCMs to better simulate the key features 
of the observed climatic episodes. In other words, if El Nino/La Nina anomaly is all one can 
hope to simulate in advance, then the most one can expect is to simulate the circulation features 
related to it. Also, one can benefit from a similar analysis of other forcing datasets such as soil 
moisture, snow/ice cover, and vegetation. 

The NDVI anomaly during the spring and summer of 1988 in the ISLSCP Initiative I data 
(Meeson et al., 1995) identifies the drought regions (Fig. 3b) experiencing a reduction. In this 
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way, the drought of 1988 provides a remarkable resemblance with these modes. This would 
suggest that the SST anomaly of 1988 had a causal influence on the production of the Great 
Plains drought. For simulating the influence of these interactions, observed SST anomalies can 
be provided to the model. If SSTs were important, a model of some reasonable credibility must 
simulate the drought. The other forcing parameter implicated in some other studies (Karl et al., 
1993) is the reduced winter snow cover in the Northern Rocky Mountains leading to low soil 
moisture and reduced Mississippi river flow with reduced irrigation and evapotranspiration 
accompanying it. One can hope to capture some of these initialization(s) through soil moisture 
and snow cover analysis produced under the Global Soil Wetness Project (GSWP, Dirmeyer et 
al., 1999), while the influence of the vegetation drought feedback can be assessed by comparing 
drought versus non-drought effects using observed vegetation data. In this way, the GEOS 
model can better capture the influence of realistic soil moisture. Atlas et al. (1993) used the GLA 
GCM to show that prescribing observed tropical SST anomalies and estimated Great Plains soil 
moisture anomalies greatly improved the 1988 U.S. drought simulation. In this paper, we also 
draw upon wind initialization using the analysis of observations produced with the GEOS 1 
version of the DAO-DAS (Schubert et al., 1993). Presumably, the recent model improvements in 
the land hydrology/snow and precipitation processes can be expected to better simulate the 1988 
drought. In addition, use of the analyzed data to better initialize the model can suggest potential 
benefits of observations to improve simulation skills. These tools and datasets provided the 
ultimate motivation for this attempt to simulate the mid-western North American drought of 
summer of 1988. 

Section 2 below describes the model used in this study. The design of the experiment is detailed 
in Section 3. Simulation results and analysis are presented in Section 4, and discussion and 
conclusions are found in Section 5. 
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2. Model Description 


The version of the GEOS GCM employed in this study had a 2 deg. (lat) X 2.5 deg. (long.) X 20- 
sigma layer resolution. The three key modules of the model are hydrodynamics, atmospheric 
physics including clouds and radiation, and Earth-Atmosphere interactions including air-sea 
interaction, biosphere, and hydrology. The hydrodynamics are on a C-grid (Takacs et al., 1994) 
with sigma layers in the vertical. This hydrodynamics has appropriate filters to eliminate 2-Ax 
modes of the dynamical atmosphere and topography (that would generate them) and the pole- 
problems. The recent developmental history of the model includes some major refinements and 
upgrades to its physical processes such as radiation and new biospheric and boundary layer 
parameterizations, as well as substantially higher horizontal and vertical resolution than used 
here. Other key features of the GEOS GCM are: i) the ability to perform coordinate translation 
and rotation with a proviso for relocating the mathematical poles to any arbitrary location (not 
used in this investigation); and ii) inclusion of a gravity-wave drag parameterization due to Zhou 
et al. (1996). Its land surface model is the so-called HY-SSiB (SSiB upgraded with hydrology 
and snow physics, Sud and Mocko, 1999 and Mocko and Sud, 2001). The convective 
parameterization of the GCM is the Microphysics of clouds with Relaxed Arakawa-Schubert 
Scheme (McRAS, Sud and Walker, 1999a&b). These packages were summarized in a recent 
paper by Sud et al. (2002). The cloud-ice fraction is diagnosed as a linear function of 
temperature - it is zero at 253. 15K and grows to unity at 233. 15K. When both ice and water 
clouds coexist, the optical thickness of the mixture is the sum of the mass fraction weighted 
optical thickness of both cloud-species. The boundary-layer scheme for turbulent transport is by 
Helfand and Lebraga (1988). The radiation package of McRAS is due to Chou and Suarez 
(1994) with a provision for handling prognostic clouds and incloud water and ice fractions (Chou 
et al., 1998 & 1999). The radiation is not too different from that of the original version of the 
GEOS GCM, except for a revised calculation for the optical thickness of clouds for short and 
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longwave radiation. For a more detailed description of different modules and parameterizations, 
the reader may refer to the original papers given as references. 

In addition to running the model with interactive physics and full dynamical responses of all the 
prognostic variables, the wind and/or soil moisture analysis data was inserted into the model for 
several cases. The wind data was taken from the DAO analysis produced at the same resolution 
and virtually using the same model. Additionally, analyzed soil moisture data, which was 
produced in a GSWP-like manner using the offline HY-SSiB model at the same 2.0 deg. x 2.5 
deg. The soil moisture data is available at three levels: surface (diurnal) layer, root zone 
(seasonal), and deep (recharge) level. The insertion of analyzed data is done in a straight-forward 
way using a direct insertion approach, i.e., simply replacing the simulated fields with the 
analyzed one at the appropriate time interval at which the analyzed data were available. The final 
two cases in this study used ISLSCP vegetation parameter data, in place of the GCM’s vegetation 
climatology. 


3. Design of the Experiment 


It is well known that most GCMs have been unsuccessful in simulating the 1988-summer (June- 
July- August, JJA) drought over North America. According to Lau’s hypothesis and data analysis 
(2001), it is partly because most scientists emphasize only the tropical sea-surface temperature 
anomalies; they often prescribe tropical anomalies as Nino 1, 2, and 3 to influence the climatic 
response in model simulations; they often ignore the influence of ocean-atmosphere interaction at 
higher latitudes as small and inconsequential. This is not an issue in this study because we 
prescribe the observed SST everywhere. In all our simulation experiments, the best estimates of 
the observed SSTs were used. Moreover, we specifically designed our simulation experiments to 
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differentially discern the influence of local, internal dynamical, and large-scale external forcings 
on the model-simulated circulation and rainfall. We have conducted seven sets of simulations as 
described below. 

1) Case 1 (Cl): a free-running model integration with prescribed SSTs, while everything else 
was fully interactive and prognostic; 

2) Case 2 (C2): model integrations ingesting analyzed winds in which simulated wind fields 
were replaced with those of DAO-analyzed winds at 6 hr intervals at all grid points; 

3) Case 3 (C3): model integrations ingesting analyzed soil moisture in which the simulated 
soil moisture fields were replaced by the analyzed soil moisture once a day at all grid 
points; 

4) Case 4 (C4): model integrations which ingested both analyzed winds and soil moisture 
following the practice in C2 and C3; 

5) Case 5 (C5): model integrations similar to C4, but with analyzed fields ingested only 
outside the limited area region shown in Figure 5; the figure also shows the La Nina cool 
anomaly over the tropical Pacific in 1988, as well as the warmer surface temperatures over 
North America in the GEOS 1 reanalysis. 

6) Case 6 (C6): model integrations in which C5 was modified with additional insertion of 
observed vegetation parameters (vegetation cover fraction, greenness, leaf area index, 
surface albedo) from the ISLSCP data (as opposed to the GCM’s climatological vegetation 
data) within the region of study; 

7) Case 7 (C7): model integrations in which C6 was updated with insertions of the opposite 
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year’s vegetation data, i.e., 1988 LAI was used for the 1987 simulation and vice versa. 


Each case contained an ensemble of four simulations that started from four consecutive days of 
December 30, 31, 1986, January 1, 2, 00Z of 1987; and December 30, 31, 1987, January 1, 2, 
00Z of 1988. Each simulation was examined from 1 June to 31 August periods of 1987 and 
1988, respectively. In highly constrained simulations, such as C2 in which simulated winds were 
replaced with the analyzed, the intra-ensemble variability was small as expected and as seen in 
the analysis of model output. For each simulation, the initial conditions of the atmosphere were 
interpolated from ECMWF analysis, whereas soil moisture and snow cover were taken from the 
GSWP-style offline HY-SSiB analysis (Sud and Mocko, 1999). In all cases, the SST (prescribed 
as monthly data) was interpolated to produce a slowly-varying daily SST using a linear 
interpolation. Therefore, Cl simulations really represent the model’s response to SST anomalies 
(warm episode of 1987 and cold episode of 1988), plus some influence of the initial soil moisture 
and snow cover prescribed at the beginning of the year. If the model were a perfect simulator of 
the Earth-Atmosphere system, one would expect the soil moisture to evolve realistically. 
However, this has not been achieved in any simulation thus far and we believe it relates to the 
unpredictability of weather on such long timescales, which affects both rainfall episodes and soil 
moisture. As a result, we did not expect the current simulation to significantly improve upon the 
model’s inability to simulate the correct land surface boundary forcings for the drought of 1988. 
Indeed, we were pleasantly surprised to find that the model did pick up some large-scale features 
(discussed in Section 4) of the circulation that reflect the existence of the drought of 1988 as seen 
in the 1988 minus 1987 JJA differences. Case C2 constraints the moisture transport but not the 
convergence, which is largely determined by the heating fields generated by the model’s physics. 
Since the soil moisture is an important forcing that is crucially affected by the precipitation and is 
the first feedback that shows large biases in response to erroneous precipitation, we decided to 
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provide the soil moisture from the HY-SSiB integration in Case C3. Thus, Case C3 generated a 
set of simulations for 1987 and 1988 in which everything was same as in Case Cl, except that 
the soil moisture was updated on a daily basis. Since these two insertions had a beneficial 
improvement, we replaced both winds and soil moisture with the analysis data in Case C4. We 
expect that if slowly-varying boundary forcings have some value, this case, with all the correct 
forcings invoked, would produce a better forecast than each of the other three: Cl, C2 and C3. 
Subsequently, we ran Case C5 to examine the influence of replacing soil moisture and wind 
outside the limited area region, and using ISLSCP vegetation data from the correct and wrong 
years in Cases C6 and Cl, respectively to find out how the biosphere feedbacks on the hydrologic 
cycle. These simulations helped us to discern: i) the factors that influenced the drought of the 
summer of 1988; ii) the influence of wind and soil moisture biases; iii) the influence of the wind 
and soil moisture biases that convey into the region from outside; and iv) the advantage of using 
observed as opposed to the climatological biosphere. We shall describe the results in the next 
section. 


4. Results 

We will describe each of the seven simulations while comparing them with one another and 
evaluating them vis-a-vis the analyzed data, so-called best estimate of observations. All of the 
analysis in Figures 7 through 12 will be shown for time-averaged June-July-August (JJA) fields. 
Since there were four cases in each set of ensemble-simulations, all results are presented as 
ensemble averages unless it is specifically stated to be otherwise. 

4.1 Case Cl: Control Simulation 

Based on the 2.0 deg. X 2.5 deg. horizontal resolution of the GEOS GCM employed for the 
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study, we would expect to simulate only the synoptic scale character of circulation changes that 
are forced by the observed SST and evolving soil moisture, snow cover and land hydrology 
anomalies. Figure 6a shows the 200 hPa streamfunction of 1988 minus 1987 for the JJA period 
simulated by the GCM (top panel) vis-a-vis the same fields from analysis of observations 
produced by DAO Data Assimilation System (DAS) (bottom panel). Positive (negative) 
differences in stream functions over the northern mid-latitudes (tropics) are evident in both plots. 
The model simulates a positive streamfunction anomaly over the drought region of North 
America (top panel) that has some synoptic scale resemblance with analyzed data for the same 
period. The similarity of these large-scale patterns indicates that the model has some skill at 
those scales. A similar examination of streamfunction differences between simulated and 
analyzed data at the 500 hPa level (Fig. 6b) again shows some resemblance between them. These 
two figures suggest that the model should show some skill in prediction of the very large-scales 
in Cl. 

Let us now examine the North American circulation and rainfall. There are large differences 
between the model and GPCP rainfall fields (Huffman et al., 1997) as shown in Fig. 7a. The 
model simulates a drought in JJA of 1988 in the mid-western to eastern United States with 1-2 
mm/day reduction in rainfall while the GPCP analysis has a less widespread reduction. On the 
other hand, the observations have many smaller scale details while the GCM does not (and at the 
resolution used, it can not be expected to). In comparing the surface to 800 hPa wind fields and 
sea-level pressure (Fig. 7b), one finds that the circulation and SLP shown in the form of 1988 
minus 1987 differences reveal that the control simulation is unable to simulate the details of the 
near-surface circulation and divergences. Consequently, all precipitation pattern and circulation 
anomalies are quite different between the model and the analysis of observations. From this, one 
infers that whatever drought is simulated in the 1988-1987 rainfall fields, the discernible 
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character of the 1988 minus 1987 circulation does not accompany it. The only influence is the 
large-scale control exerted through the observed SSTs, evidenced in the 200 and 500 hPa 
anomaly patterns. The model does simulate a drought in 1988 that also affects the 
evapotranspiration (Fig. 7c) and surface temperature (Fig. 7d). Both of these fields indicate that 
the GCM is drier and warmer than observed; most likely, it is a consequence of a positive 
feedback between soil moisture and surface temperature. Excessive evapotranspiration is likely 
to cause decreasing soil moisture and higher temperatures, particularly in the summer. However, 
it would appear that the model does simulate some sort of a large-scale drought in JJA of 1988 
with respect to 1987. Such a forecast could be useful, but it has very little scientific value for a 
researcher because of the missing details and biases in large-scale anomalies. 

4.2 Case C2: Wind-Updated From Analysis 

We next examine the influence of local, synoptic scale and collateral errors on the JJA 
circulation and rainfall. This is done in the simulation experiments C2 through Cl. Case C2 
represents a simulation in which the simulated winds are replaced by the observed at 6-h 
intervals. Any horizontal wind field has two components: one is divergent and the other is 
rotational. The updated winds would alter the rotational part, but the divergent part will be 
determined by the model’s diabatic heating fields, particularly the temperature (which is not 
updated) and its effect on associated divergence. As far as the rotational part is concerned, these 
winds do not change much, but help to transport heat and moisture as well as alter the pressure 
gradients to establish the geostrophic (vorticity-pressure) relationship. With observed winds, the 
model must capture the observed transports while it can modulate its divergences in response to 
diabatic heating fields produced by the model’s physical interactions. Over a short time-period, 
dynamics generally overwhelms the physics; therefore, by inputting observed winds, the 
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dynamics gets constrained everywhere, whereas thermodynamics has all the freedom to operate 
according to the model’s physical parameterizations. 

Naturally, the C2 simulated 1988 minus 1987 precipitation fields (Fig. 8a) and sea-level pressure 
(Fig. 8b) (surface to 800 hPa motion fields as prescribed) are in much better agreement with 
observations. There was virtually no difference in any fields including precipitation among the 
members of the ensembles. The evapotranspiration patterns (Fig. 8c) also had a strong 
resemblance among the members. Any differences can be either due to soil moisture anomalies, 
or radiative forcing or any remaining SSiB deficiencies in the simulation of the entire biosphere 
land-hydrology complex. Comparison of differences over the ocean with offline Global Soil 
Wetness Project (GSWP) data is not meaningful because GSWP does not apply over the oceans. 
The accompanying surface temperature differences (Fig. 8d) are consistent with the 
evapotranspiration differences, with somewhat better resemblance to GSWP and GEOS 1 
reanalysis fields (Figs. 7c & 7d). Comparison of simulated soil moisture to the GSWP soil 
moisture for this case found that the soil moisture was generally drier for C2 than for Cl. 
Overall, we note that wind errors in the chosen region are responsible for most of the synoptic 
scale errors. This suggests that one can benefit substantially by having a forecast system in 
which winds are properly initialized or better simulated. This is in agreement with a number of 
previous studies (e.g., Atlas et al., 2001) in which importance of winds is considered vital to the 
accuracy of forecasts. Any forcing that modulates winds will affect the forecasts. 

4.3 Case C3: Soil-Moisture Updates from GSWP Analysis 

Instead of winds, here we used the GSWP-analyzed soil moisture, which was performed by 
replacement of simulated soil moisture with the offline analysis (add reference) at each grid-point 
at 00Z (once a day). A cursory glance at the precipitation, sea-level pressure and surface to 800 
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hPa winds, evapotranspiration and surface temperature fields (Figs. 9a-d) shows that it did not 
have much of an impact beyond making the evapotranspiration anomalies slightly better than 
Case Cl. The model’s patterns represent a combination of soil moisture and net radiative forcing 
and circulation. At least over the drought regions, reduced soil moisture has produced warmer 
temperatures even though the synoptic scale surface to 800 hPa circulation remained largely 
unaffected. 

4.4 Case 4: Both Wind and Soil Moisture Updates 

In this simulation, both winds and soil moisture were prescribed from analysis of observations. 
Precipitation differences (Fig. 10a) have a much closer resemblance with GPCP rainfall 
differences (Fig 7a). As expected, most of the benefits were derived from input of observed 
winds. Soil moisture benefits were much smaller. The SLP changes are similar to those of C2 
and the influence of analyzed soil moisture in not discernible in the JJA average. However, in 
the case of evapotranspiration anomalies (Fig. 10c), the combined simulation is similar to the soil 
moisture anomaly simulation C3 (Fig. 9c). This shows that for evapotranspiration, which is a 
major component of surface energy partitioned into different surface fluxes, soil moisture must 
be more important than the wind. In that way, the current result makes good intuitive sense. A 
clear difference is in the surface temperature anomaly patterns. They are similar to the surface 
temperature patterns of C3 (Fig. 9d), but the intensity of differences is much reduced and the 
values are in better agreement with the analysis. In this way, the use of observed winds helps 
transport the air mass and its associated temperature more realistically as compared to the soil 
moisture update only, Case C3. 

4.5 Case C5: Winds and Soil Moisture Updates Outside LBOX 
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Case 5 represents simulations with both soil moisture and wind fields updated, but only outside 
of the region of 23-6 IN and 66-129W (hereafter, the LBOX region). It is a simulation in which 
the influence of all external-forcing errors in the two chosen fields is removed by having the 
influx of energy and water vapor outside the chosen region to get updated with analyzed winds 
and soil moisture. This in essence is equivalent to running a regional model with the best 
available input of humidity and winds. Therefore, Case C5 must be compared to C4 and Cl. 
Case C5 has large similarities with Case C4 and in some respects it is even a better simulation 
than all the others. Figure 11a shows that the precipitation field gives a better simulation of the 
magnitude of the observed drought extending from Canada to the north of the Great Lakes 
through Wyoming, although the orientation of the drought west of 90W is somewhat poorer. The 
magnitudes are somewhat smaller than that of C4, which is an improvement. The primary mid- 
western drought region has synoptic-scale character, but it is equally well/poorly simulated in 
both C4 and C5, i.e., the deficiencies and strengths of both similar. The only difference is that 
magnitudes of Case C5 are smaller and are in better agreement with data. The changes in other 
areas are unremarkable. In C2 and C4, the circulation was prescribed, but in C5 it is only 
prescribed outside the box; consequently, its SLPs (Fig. lib) were degraded as compared to C4 
or C2. This implies that model-introduced errors inside the dynamically free box make SLPs 
drift away from observations (as expected). The case with prescribed winds and soil moisture 
(C4) had very little inter-ensemble variability, whereas Case C5 has much more (although not 
nearly as much as Cl). This drift is related to the model’s freedom to evolve its own circulation 
and hydrologic processes in the region. Figure 11c shows evapotranspiration anomalies. GSWP 
evapotranspiration anomalies (Fig. 7c) are much smaller than that of Cases C4 and C5. In the 
higher latitudes, where there is enough soil moisture, the 1987 minus 1988 JJA 
evapotranspiration anomalies are not so large, but in the mid- western drought regions they follow 
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the observed precipitation forced initial soil moisture for Case C4 and simulated precipitation for 
Case C5. The surface temperature anomalies (Fig. lid) are dependent upon winds, cloudiness 
(affecting solar radiation reaching the surface as well as net outgoing longwave radiation), and 
evapotranspiration. Since winds are prescribed, the only remaining degrees of freedoms are 
cloudiness and soil moisture, which produce the observed effects. 

4.6 Cases C6 and C7: Same as Case C5 with observed vegetation inside LB OX 

In view of a number of sensitivity studies highlighting the importance of biosphere-atmosphere 
interaction, here we examine how useful is GEOS model’s sensitivity in simulating the drought 
circulation. Figure 12a shows that with actual vegetation data, the drought in the mid- west shrank 
somewhat more than that of Case C5 as compared to analysis of observations. In addition, the 
East Coast of North America was wetter than that of C4, which itself was wetter than rainfall in 
GPCP data (Fig. 7a). In this respect, the observed vegetation did not help. In the SLP fields 
(Figure 12b), the differences between C5 and C6 are unremarkable. The evapotranspiration 
anomalies in Figure 12c mimic rainfall anomalies, suggesting that if one simulates large errors in 
the rainfall, the evapotranspiration through soil moisture feedback will change correspondingly 
regardless of vegetation parameters. In this simulation many other parameters that are associated 
with the modified vegetation could not be realistically altered. However, since the parameters 
modified are considered to be the dominant modulators of evapotranspiration, this should not 
make much difference. It would be expected that the drought vegetation parameters of C6, 
which are less than in C5, will cause less evapotranspiration; however, even this does not happen 
because the biospheric feedback interactions are so complex that changing the parameters did not 
affect the time-mean rainfall correctly to make much difference to the simulation. On the other 
hand, rainier East coast produced higher evapotranspiration and cooler surface temperatures. 
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Case Cl (with opposite year’s vegetation data) was very similar to C6 as far as JJA 1988 minus 
1987 rainfall, SLP and winds from surface to 800 hPa, evapotranspiration and surface 
temperatures (Figs. 13 a-d) showing that in our study, there were no discernible differences in the 
simulation as a consequence of the correct or incorrect vegetation data for 1987 and 1988. It 
should be evident from the comparison between C5 and C6, but the Cl simulation helps to 
reinforce the result, mainly as a result of the ISLSCP vegetation data for 1987 and 1988 being 
closer in agreement to each other than either was to the vegetation data used in the GEOS GCM’s 
climatology. Dirmeyer (2000) with the COLA GCM, which has essentially the same SSiB, has 
shown that correct soil moisture helps in 1987 and 1988 switched simulations; indeed, it does so 
in our case too (not shown), but the correct vegetation effect is really small as compared to the 
other parameters. 

4.7 Analysis of Individual Cases in the Ensemble 

The 1988 minus 1987 JJA precipitation for the four individual ensemble members for Case Cl 
are shown in Figure 14a. The right two panels of this figure show the mean of the four members 
(top) and the GPCP analysis (bottom). For simulations Cl-a and Cl-d, the drought is simulated 
mainly over the Midwest and eastern U.S. However, the other two ensemble members simulate a 
widespread drought also over the Great Plains across the northern Rockies. This figure 
demonstrates both the ability of the GCM to predict a drought six to eight months in advance, as 
well as the uncertainty in predicting its location as a result of differing initial atmospheric states. 

For Case C4, Figure 14b shows how strongly constrained the precipitation anomaly becomes by 
just replacing the simulated wind and soil moisture fields with the analyzed wind and soil 
moisture fields. Very minor differences can be noted, but the general pattern of a drought in the 
upper Midwest to Great Plains and around the Great Lakes into eastern Canada is virtually 
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identical among all four ensemble members as well as the GPCP panel. For Case C5 (Fig. 14c), 
with winds and soil moisture only replaced outside the region, greater inter-ensemble differences 
are found, but not nearly as large as in Case Cl. 

The individual ensemble members were also analyzed with a cyclone tracking routine of Terry 
and Atlas (1996) for JJA for Cases Cl, C4, and the reanalysis (not shown). Results from the 
reanalysis showed no significant difference in the location or frequency of cyclones between 
1987 and 1988 during this period. Furthermore, almost no cyclone activity was identified in the 
northern Great Plains during both years. Similar results were found for Case C4, with all four 
ensemble members having very similar cyclone tracks, as a result of the replaced winds. In the 
control Case Cl, several cyclones among all ensemble members were noted in this region during 
both years, as well as considerable scatter between individual members. 

4.8 Behavior of Regional Averages 

Monthly- averaged precipitation plots for Cases C1-C6 for the Great Plains region of 30-50N and 
85-100W for 1988 are shown in Figure 15a. The solid line represents the GPCP data and the 
thick dashed line represents the ensemble average for that case. The four thin dashed lines are 
the precipitation for each individual ensemble member. For Case Cl, large differences between 
ensemble members were produced, but the average shows that on the whole, the GCM failed to 
both simulate the strong spring and early summer drying, as well as the late summer return of the 
precipitation in this sub-region. The GCM was too wet in the early summer largely due to a 
poorly simulated circulation, while it was too dry in the late summer from a positive feedback of 
progressively lower soil moisture in this region. For Cases C2 and C4, the strong inter-ensemble 
member similarity is again noted, with the addition of soil moisture data in C4 helping to 
simulate better late summer precipitation. The soil moisture also had a positive effect on 
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simulated precipitation in C3. Case C5 also shows the problem with late summer drying 
affecting the precipitation, with the ensemble members being much more similar to each other. 
Adding the correct vegetation data in Case C6 only somewhat reduced the anomalous high early 
spring precipitation and low late summer precipitation also shown in C5. Figure 15b shows the 
same data for 1987 with the spring and early summer precipitation drought not found in the 
observations. Results from the GCM are generally similar to 1988, however, with the addition of 
soil moisture data helping the simulation of late summer precipitation, and the box region tightly 
constraining the simulations, but still not simulating the situation correctly. 

Daily-averaged soil moisture data for the Illinois region bounded by 38-4 IN and 88-92W for 
1988 are shown in Figure 16a. Here the solid line represents the soil wetness from the offline 
HY-SSiB analysis forced with the ISLSCP Initiative I data (Sud and Mocko, 1999), which had 
been previously shown to well reproduce observations of soil moisture in this region. Again, a 
large scatter is evident amongst ensemble members in the control Case Cl. When the winds are 
replaced in Case C2, the stronger than observed late summer drying of the soil is seen. The 
simulated soil moisture is very close to the observed in Cases C3 and C4 as a result of replacing 
soil moisture data daily, but the effect of free-running winds is noted with the numerous soil 
moisture spikes in C3 before the 00Z replacement. Cases C5 and C6 agree with previous results, 
with a moderate amount of scatter and error. In Figure 16b, the soil moisture feedback error is 
further highlighted. In Case C2, the late summer soil moisture is much drier than observed, as 
are Cl, C5, and C6. 

5. Discussion and Conclusion 

As discussed in the Introduction, the drought of 1988 had not been simulated satisfactorily by 
any general circulation model. Its physical description, together with the available data analysis 
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to-date, shows that its primary cause was the SST anomalies of 1988. In fact, 40 years of data 
analysis has isolated the very same regions of North America for the occurrence of droughts in 
response to La Nina SST anomalies. Therefore, the JJA 1988 was a very typical drought and one 
would expect a GCM to be able to simulate it. 

When the SST anomalies were prescribed, the GEOS GCM did produce some of the circulation 
features of a drought over North America, but these features could be identified only on the 
planetary scales. The 1988 minus 1987 precipitation fields show that the GCM was successful in 
simulating a drought for the JJA period with reduced precipitation, but the accompanying 
circulation was so poor that one would probably infer that the GCM simulated the drought for the 
wrong reason. To isolate the causes for the above behavior, winds and soil moisture were 
prescribed from analyses of observations as continuous updates to the simulation. Other fields 
such as temperature, humidity, and/or surface pressure could not be used without invoking full 
data assimilation. A simple reason is that any large mismatch of incoming and simulated data 
must be dealt with special data inconsistency-resolving techniques. If the offline model’s snow 
cover, for example, does not match that of the new simulation, the huge temperature difference 
between the original skin temperature and inserted air temperature after the direct insertion of 
data would cause boundary layer fluxes to respond in an outrageous fashion and force the 
boundary-layer code to become unstable. Wind fields and soil moisture fields are not as 
disruptive and could be easily assimilated by direct replacement either individually or together as 
was done for Cases C2 through C4. 

Our results show that a large number of simulation biases emanate from wind biases that are 
carried into North American region from surroundings areas. When using analyzed winds, the 
rotational part of the circulation is prescribed from analysis, whereas the divergent part is 
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modulated by the model’s physics. With winds prescribed at 6-h intervals, the inter-ensemble 
variability of the simulations virtually vanishes. In addition, assimilated winds produce a much 
better simulation of both the precipitation and low-level circulation at the model s resolution. 
Inclusion of soil moisture also helps to ameliorate the strong feedback, perhaps stronger than that 
of the real world, between soil moisture and precipitation that produced too strong precipitation 
anomalies in the control case. Whatever remaining differences between the observed and 
simulated precipitation and surface temperatures exist, they are caused by errors in the model s 
physics which includes the cloud radiation interaction, the precipitation process itself, and the 
land-atmosphere interaction. The simulations showed the structure of surface temperature and 
precipitation errors in response to winds alone, soil moisture alone, and both. For the case of 
prescribed winds, the surface temperature anomalies have one persistent pattern, whereas for the 
soil moisture it is another. In the combined case, the two patterns merge and also give somewhat 
more realistic evapotranspiration and precipitation patterns. 

The cyclone track analysis did not show a meaningful difference between 1987 and 1988 for JJA, 
in either the GCM or the reanalysis. These cyclones tend to be weaker and less frequent in the 
summer months, and both the observations and model show a strong precipitation deficit in 1988 
despite little change in cyclones. Thus, the cyclones likely produced only a small amount of the 
precipitation in the Great Plains, which is in agreement with Fritsch et al. (1986) who showed 
that mesoscale convective systems (MCSs) account for 30 to 70 percent of the summertime 
precipitation. In the current configuration, the GCM is unable to resolve MCSs, and thus we 
were unable to examine the strength and frequency of MCSs in this region for both years. 

Case C5 with the LBOX region without updating soil moisture and winds inside the region, 
while outside the region winds and soil moisture were updated as in Case C4, showed the 
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following. Case C5 is virtually of the same quality as Case C4, and shows that most of the local 
simulation errors originate outside the LBOX region. This is a not-unexpected result because if 
weather and climate have global connectivity, then any chaotic component of weather will 
propagate into a region such as LBOX from outside; however, such a large magnitude of this 
connectivity, even on a seasonal scale, is a new result. Even though we do not know about the 
robustness of this finding for other models, one naturally expects it would not be too different for 
other state-of-the-art climate models. The new result also reaffirms how and why regional 
models run with prescribed lateral transports from data are able to do a more realistic job of 
simulating a specific phenomenon while a global model continues to have problems. However, 
in the long run, only a free-running GCM will enable scientists to simulate climate for the future. 
In that spirit, this research is not an end in itself, but helps to provide guidance as to the 
weaknesses of the GCM. One naturally wonders if weather is not deterministically predictable 
beyond 5-10 days, will its time-mean (climate) also contain a significant component of 
unpredictably over the three months (June, July, and August) being analyzed. The question boils 
down to finding out if the model’s biases, which also contribute to the lack of predictability, are 
so large as to limit the value of its predictions. On the other hand, the model is not able to 
capture any discernible improvement in simulations with the correct versus opposite vegetation 
parameters. Clearly, vegetation and soil moisture go together, but if the dominant role is of soil 
moisture, then its biases would mask any plausible advantage of using the correct vegetation. In 
a similar study by Dirmeyer (2000), when only the correct versus the opposite soil moisture was 
provided, it revealed a discernible improvement in response to correct soil moisture. The new 
results show that the soil moisture governs the outcome of land atmosphere feedback interaction 
much more than the vegetation parameters. 

The model’s biases in the prescribed soil moisture simulations as well as in the prescribed wind 


25 


simulations are quite persistent. The authors believe that model biases significantly influence the 
prediction errors. These biases are so dominant that through wind-errors, they change the 
transports of heat and moisture into the LBOX region. When winds are prescribed from analysis 
of observations outside the chosen LBOX region, the model produces a much better JJA drought 
as compared to Control, and it remains almost as good as the one in which winds and soil 
moisture were prescribed everywhere. This shows that biases in circulation and advective 
transports propagate and strongly contribute to the simulation biases. SCM simulation research 
(Ghan et. al., 1999 and Xie et al., 2002) shows significant model-physics-dependent biases 
among the participating models even over a single grid-cell; therefore, the authors conclude that 
the only meaningful way to improve these GCMs is to first reduce their biases at the grid-cell 
level. This would require improvements in cloud physics, cloud radiation interactions, boundary- 
layer processes, as well as the rest of the atmospheric column physics. Without such a concerted 
effort in model improvement, simulating climatic events will continue to be a hit-or-miss 
prediction. Consequently, its global change inferences will continue to be unreliable. 
Furthermore, even when the model realistically forecasts a climate event, scientists will ponder 
about the right/wrong reason for the forecast, and that in turn will haunt modelers attempting to 
simulate global change scenarios. 

In the series of simulation attempts reported above, the model is able to predict the very large 
scale circulation changes as seen in the 200 hPa streamfunction differences of 1988 minus 1987 
JJA circulation somewhat reasonably over North America, which leads to a somewhat misplaced 
Midwestern drought in 1988 minus 1987 rainfall. However, the rest of the circulation is not well 
reproduced. Presumably, soil moisture, which has a much longer time scale (but can be affected 
by a single weather event, whose course is largely unpredictable), is modulating the ensuing 
circulation. Even when the soil moisture is prescribed from GSWP analysis of 
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hydrometeorological data from analysis of observations for 1987 through 1988, the 
evapotranspiration errors remain large. This implies that the net radiation at the surface and 
vertical temperature and humidity structures that are governed by thermodynamics and vertical 
column adjustment physics of the model are contributing to the biases. In fact, if cloud 
distribution or cloud radiative feedback are erroneous, net radiation at the surface would be 
affected and that will influence the evapotranspiration and Bowen ratio. It appears there are 
significant modeling errors associated with not being able to simulate the drought well, even 
when winds and soil moistures are prescribed everywhere. 
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Figure Captions 


Figure 1: Meteorological aspects of the North American drought of 1988. Figure adapted 
from NCEP/NWS. 

Figure 2: 1988-1987 June/July/August (JJA) Sea Level Pressure (in mb) and surface to 800 
hPa averaged wind vectors (in m sec -1 ) from the GEOS 1 Reanalysis. 

Figure 3: NDVI anomaly for August 1988 over North America from: a) data gathered 
wdth NOAA’s Advanced Very High Resolution Radiometer (AVHRR); and b) data from the 
ISLSCP Initiative I Dataset at 1 degree global resolution. 

Figure 4: Typical North American pattern during the summer months associated with a 
cold event in the eastern Pacific (La Nina). Figure adapted from COAPS/Florida State 
University. 

Figure 5: 1988-1987 June/July/August (JJA) Surface Air Temperature (in K) from the 
GEOS 1 Reanalysis. The region inside the rectangle over North America from 23-61 North 
and 129-66 West is allowed to freely vary as described in the LBOX experiment. 

Figure 6a: 1988-1987 June/July/August (JJA) 200 hPa streamfunction (in 10 10 kg sec -1 ) 
from: top) Case Cl CTRL; and bottom) GEOS 1 Reanalysis. 

Figure 6b: Same as Figure 6a, only for 500 hPa streamfunction. 

Figure 7a: 1988-1987 June/July/August (JJA) precipitation (in mm day -1 ) from: top) Case 
Cl CTRL; and bottom) GPCP Version 2. 

Figure 7b: 1988-1987 June/July/August (JJA) Sea Level Pressure (in mb) and surface to 
800 hPa averaged wind vectors (in m sec -1 ) from: top) Case Cl CTRL; and bottom) GEOS 
1 Reanalysis. 

Figure 7c: 1988-1987 June/July/August (JJA) evaporation (in mm day -1 ) from: top) Case 
Cl CTRL; and bottom) GSWP offline. 

Figure 7d: Same as Figure 7b, only for surface air temperature (in K). 

Figure 8: Same as Figure 7, only for Case C2 WIND. 

Figure 9: Same as Figure 7, only for Case C3 SOIL. 

Figure 10: Same as Figure 7, only for Case C4 BOTH. 

Figure 11: Same as Figure 7, only for Case C5 LBOX. 
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